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ABSTRACT 

Aims. We investigate the surroundings of the hypercompact Hu region M 17 UC1 to probe the physical properties of the associated 
young stellar objects and the environment of massive star formation. 

Methods. We use diffraction-limited near-IR (VLT/NACO) and mid-IR (VLT/VISIR) images to reveal the different morphology at var¬ 
ious wavelengths. Likewise we investigate the stellar and nebular content of the region by VLT/SINFONI integral field spectroscopy 
with a resolution R ~ 1500 at H + K bands. 

Results. Five of the seven point sources in this region show L-band excess emission. Geometric match is found between the FL 
emission and near-IR polarized light in the vicinity of IRS5A, and between the diffuse mid-IR emission and near-IR polarization 
north of UC1. The FL emission is typical for dense PDRs, which are FUV pumped initially and repopulated by collisional de¬ 
excitation. The co-presence of Hei, Hi, and H 2 lines in most region argues against an edge-on configuration of the M 17 SW PDR, 
but is in favor of a moderately inclined geometry with respect to the line of sight. The spectral types of IRS5A and B273A are B3-B7 
V/III and G4-G5 III, respectively. The observed infrared luminosity L lR in the range 1-20 pm is derived for three objects; we obtain 
2.0 x 10 3 L 0 for IRS5A. 13 L 0 for IRS5C, and 10 L e for B273A. 

Conclusions. IRS5 might be a young quadruple system. Its primary star IRS5A is confirmed to be a high-mass protostellar object 
(~ 9 M 0 , ~ 1 x 10 s yrs); it might have terminated accretion due to the feedback from the stellar activities (radiation pressure, outflow) 
and the expanding Hii region of M 17. UC1 might also have terminated accretion because of the expanding hypercompact Hu region 
ionized by itself. The disk clearing process of the low-mass YSOs in this region might be accelerated by the expanding Hu region. 
The outflows driven by UC1 are running in south-north with its northeastern side suppressed by the expanding ionization front of 
M 17; the blue-shifted outflow lobe of IRS5A is seen in two types of tracers along the same line of sight in the form of H 2 emission 
filament and mid-emission. The H 2 line ratios probe the properties of M 17 SW PDR, which is confirmed to have a clumpy structure 
with two temperature distributions: warm, dense molecular clumps with n H > 10 5 cnr 3 and T ~ 575 K and cooler atomic gas with 
n H ~ 3.7 x 10 3 - 1.5 x 10 4 cm- 3 and T ~ 50 - 200 K. 

Key words, stars: early-type -stars: formation -stars: individual: M 17 UC1 -ISM: individual objects: M 17 SW -ISM: molecules: 
photon-dominated region (PDR) 


1. Introduction 


Massive stars (> 8 M 0 ) affect their surroundings by ionizing ra¬ 
diation and strong stellar winds throughout their life, as well 
as by metal enrichment in their final fate as Supernovae. How¬ 
ever, the formation process of massive stars is still a hot open 
debate. Both theoretical prediction and observational evidence 
suggest that massive stars most likely form via accretion of 
material quite similar to low-mass stars, but the involved pro¬ 
cesses may not be only simply scaled up (IZinnecker & Yorkel 
120071 and references therein). The major difference is the ener¬ 
getic feedback (radiation pressure, stellar winds and outflows) 
which rapidly dissipates the circumstellar envelope and conse¬ 
quently limits the mass growth of the central young stellar ob¬ 
ject (YSO). In addition, for high-mass YSOs the earliest evolu- 


* Based on observations by the European Southern Observatory Very 
Large Telescope on Cerro Paranal, Chile (ESO program IDs: 281.C- 
5027(A), 281.C-5051(A,B)) 


tionary phase (protostellar phase, before reaching the zero-age 
main sequence, hereafter ZAMS) is very short, as characterized 
by Kelvin-Helmholtz timescales of less than 10 5 yrs. Thus, mas¬ 
sive stars evolve fast, even during their accretion phase. It is 
thought that the time spent in the main accretion phase might be 
significantly less than the main-sequence lifetime, presumably 
of the order of a few dynamical times of the star-forming molec¬ 
ular core. Therefore, to accumulate sufficiently large amounts 
of material during the very short timescale, accretion rates of 
high-mass YSOs must be much higher than those of low-mass 
YSOs, e.g., > lO _4 M 0 yr _1 (e.g. lHosokawa et al.l2010i hereafte r 
H+10) compared with < 10 _8 M o yr _1 (e.g. Fang et al1l2013l) . 
The strong feedback from massive protostars may prevent the 
accretion even before the arrival at the ZAMS if accretion rates 
are constantly larger than a few 10 3 M 0 yr 1 . This settles an up¬ 
per limit for the protostars around several tens of solar masses. 
In fact, accretion rates may vary strongly with the evolution of 
massive protostars. Massive accretion might continue in a non- 
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steady fashion, which potentially allows even more massive stars 
to form by mass accretion (H+10). 

The evolution of massive protostars at such high accre¬ 
tion rates is still the vague feature of massive star formation 
which determines the feedback on their environment. Numer¬ 
ical simulations show that accreting massive protostars with 
high accretion rates have large radius. For instance, at an ac¬ 
cretion rate of 1 x 10 3 M 0 yr 1 , the protostellar radius may 
exceed_100 R Q _at maximum (e.g. lYorke & Bodenheimerll20(M 
lHosokawa & Omukaill20Q9l) . Such a large radius hence leads to 
low effective temperature, and very low stellar UV luminosity, 
which might b e too low for the growth of an H n region aroun d 
the protostar (IChurchwelll2002t lKetol2007l iMurphv et al.l20 1 Oh . 
High-mass protostars are considered to precede the formation of 
an H ii region. From the observational view, a lot of high-mass 
protostar candidates ha ve high infrared luminosities without 
observable H n regions dKumar & Gravel 12007b iMolinari et al.1 
(20081 lG rave & Kumarll2Q09l). 

Sect. |2] briefly introduces the massive star-forming region 
M 17, and particularly summarizes the studies of the M 17 UC- 
IRS5 region. In Sect. [3 the observations and data reduction are 
described. The results based on near- to mid-IR data are pre¬ 
sented in Sect. [4] The derived properties are discussed in Sect. [5] 
and the conclusions are presented in Sect. [6] 


2. M17 UC1 - and the southwestern 
photodissociation region 


In a distance of 1,98/(jp kpc ( Xu et a l. 201 1), M 17 is among 
the best laboratories in the Galaxy for investigating the for¬ 
mation of massive stars. In this paper, we report near- to 
mid-IR imaging and integral-field spectroscopic studies for 
the M 17 UC1-IRS5 region (see Fig. |T}. which is located 
just west of the arc-like ionization front (IF). This region is 
well known because of the hypercompact Hii (HCHIIregion 
M 17 UC1 (tSewilo et al. 2004), which is surrounded by a cir- 
cumstellar disk dNielbock et all |2007( hereafter N+07). The 
other interesting object, M 17 IRS5, is a bright IR source located 
5" southwest of M 17 UC1; in contrast to the HCHII region it 
is not detectable at 1.3 cm (iChini et al.l l2000t) . Its spectral en¬ 
ergy distribution (SED) at IR wavelengths suggests a warmer 
component with color temperature 1000 K, and a cooler com¬ 
ponent ~ 150K dNielbock et al.l 1200 ll hereafter N+01). The 
non-detection of an associated Hu region is reminiscent of 
the early protostellar phase when the protostar is huge with 
low effective temperature and low UV luminosity - meaning 
that IRS 5 might be younger than UC1. However, this expla¬ 
nation isjust one of the three plausible scenarios proposed by 
iKassis et al.l d2002l) (hereafter K+02). An H ii region with a den¬ 
sity higher than 3 x 10 5 cm 3 or a heavy dusty envelope could 
also escape detection. Previous near-IR polarization studies re¬ 
vealed infrared reflection nebulae (IRN) associated with the two 
sources (iChen et al.ll2012[ hereafter CZ+12), which might trace 
potential outflows. The third bright IR source, M 17 B273, pro¬ 
jected against the edge of the arc- shaped IF, might be a YSO too 
(near-IR excess, see Hanson et alii 1997 1 but requires a spectro¬ 
scopic classification. 

Besides characterizing the stellar content, the SINFONI 
integral-field spectroscopy can also be used to investigate 
the diffuse nebular emission of M 17 SW-one of the best- 
studied dense PDRs in the Galaxy. At far-IR to millime¬ 
ter wavelengths, studies of molecular and atomic emission 
indicate that the structure of the gas is highly clumped 


(Stutz ki et al.l 1 1 988b Stutzk i & Guestenl Il990b iMeixner et al.l 
119921: iPerez-Beaunuits et al.l 20101) . and supp orted by mag- 
netic field rather than by thermal gas pressure ((Pellegrini et al.l 
120071). Temperatures of ~ 275 K were found toward the IF 
dBro'gan & Trolandi 1200 ll) . One characteristic of PDRs is H 2 
emission originating from the collisional deexcitation of H 2 
molecules initially excited by UV photons, which js an impor- 
tant h e ating mechanism o f dense PDRs ((Sternberg & Dalgarno 
119891) , ISheffer & Wolfirel (1201 3l) reported four mid-IR pure- 
rotational H 2 lines toward M 17 SW which are consistent with H 2 
emission in high-density clumps («h > 10 5 cm~ 3 ) embedded in 
an interclump a tomic gas of density l ower by two or three order s 
of magnitude (IMe ixner et al.l 119921: IMeixner & Tielensl Il993l) . 
iNiirnberger et al.l (1 20071) reported near-IR H 2 emission associ¬ 
ated with a jet ejected by a forming high-mass protostar, which 
is just located V southeast to the M17 UC1-IRS5 region. In 
such case, H 2 emission is produced by thermal emission in shock 
fronts. 


3. Observations and data reduction 

3.1. SINFONI observations 

High spatial resolution (AO supported), middle spectral reso¬ 
lution (R ~ 1500) near-IR integral field spectroscopic data of 
the M 17 UC1-IRS5 region (see Fig. Q} were taken in Ser¬ 
vice Mode during the nights 2008-06-08/09,2008-09-26/27 and 
2008-09-27/28, using ESO’s near-IR integral field spectrograph 
SINFONI (Spectrograph for IN tegral Yield Observations in the 
Nc«r-Infr«ref/: [Gillessen et al.l2005l) mounted on the Cassegrain 
focus of the VLT Yepun at the Paranal observatory, Chile. 

Overall, we observed a sequence of 4 adjacent target 
positions (see Table |T}. In all cases, AO curvature sens¬ 
ing was performed on the reference source M17-CEN64 
(RA= 18:20:25.71, DEC = -16:11:41.7, J2000), which is lo¬ 
cated at the distance of about 1675 toward the east of IRS5 
(about 12" toward the south-east of B273). For appropriate sky 
subtraction, a nearby empty sky position was available at the dis¬ 
tance of about 70" - 90" toward the west of IRS5. 

In a trade-off between avoiding saturation of the spectra of 
the central point source and optimizing both dynamic range 
and field of view (FOV) for the diffuse circumstellar emis¬ 
sion, we set up SINFONI with its 100 mas pixel scale, cover¬ 
ing an instantaneous FOV of about 3" x3" with 64x32 pixels 
of 50 mas x 100 mas each, together with the H + K grating, 
which comprises the wavelength range from about 1.45 pm to 
2.45 pm at the spectral resolution of about 1.0 nm (dispersion of 
0.50nm/pixel). The detector integration times (DITs) were set 
to 60s (with NDIT = 1), both for on-source (science) exposures 
and off-source (sky) exposures. 

For each target position, we applied a sequence of dither off¬ 
sets on a 3 x 3 position grid with grid spacings of 075, resulting 
in an effective FOV of about 4" x 4". Because the faint circum¬ 
stellar material around IRS5 is rather widespread, we took ad¬ 
ditional exposures on a 3 x 3 position grid with grid spacings of 
270 and 272, increasing the effective FOV to about 774 x 774. In 
total, we gathered 27 exposures on IRS5, 6 exposures on the po¬ 
sition intermediate between IRS5 and UC1, 6 exposures on UC1 
and 18 exposures on B273. To allow proper sky subtraction, sci¬ 
ence exposures (with the AO loop closed, of course) were in¬ 
terleaved every 10 minutes by several (typically 5-6) exposures 
on the sky position (although void of stars and diffuse emission, 
with the AO loop open). 
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Table 1. SINFONI observation parameters 


Target 

R.A. 

(J2000) 

DEC. 

(J2000) 

Exp. Time 

(s) 

Telluric STD 

M17-IRS5 

18 20 24.60 

-16 11 39.4 

27 x60 

HIP 091126 (G2 V) 


18 20 24.71 

-16 11 37.1 

6x60 

HIP 094378 (B5V) 

M17-UC1 

18 20 24.82 

-16 11 34.9 

6x60 

HIP 094378 (B5V) 

M17-B273 

18 20 25.07 

-16 11 33.9 

18 x60 

HIP 092470 (B2IV) 


To correct for telluric features and to flux-calibrate the sci¬ 
ence exposures, each night they were immediately followed by 
observations of a suitable telluric standard star, applying the 
same strategy (auto-jitter pattern, fixed sky offset) as for the sci¬ 
ence targets and matching their airmass. The telluric standard 
stars used in the three observation nights are listed in Table [T] 

During all three nights the observing conditions were good, 
with sky transparency clear (exception: thin cirrus clouds pass¬ 
ing through during the night 2008-09-27/28) and telescope guide 
probe seeing measurements typically in the range 076 - 172. 
Correspondingly, the achieved strehl ratios were in the range 
20 - 40%, as directly measured on the spectrally collapsed data 
cubes. 

Dark frames, lamp flats, and arcs (taken with SINFONI’s 
internal Neon and Argon lamp for the purpose of wavelength 
calibration) were obtained through the SINFONI scientific cali¬ 
bration plan. All basic steps of data reduction (flat-fielding, sky 
subtraction, bad pixel and atmospheric distortion correction) as 
well as the wavelength calibration were performed with the SIN¬ 
FONI data reduction pipeline. Final merging of the fully reduced 
and wavelength calibrated individual data cubes to one mosaic 
covering all science targets was performed within the SINFONI 
pipeline, too. 

3.2. Ancillary near- to mid-IR data 

The JHKL AO imaging w as carried out in 2003 June us¬ 
ing N AOS/CONICA (NACO: iLenzen et al.ll200l iRousset et all 
l2003h on the ESO VLT at the Paranal observatory, Chile. The 
FOV is 27" x 27" with pixel resolution of 07027. The 3<x limit¬ 
ing magnitudes are J = 20.2, H = 19.7, K = 19.3, and L - 15.2. 
The photometry was carried out merely for the point sources 
within the FOV of SINFONI data (see Fig. [Q- The magnitudes of 
these point sources were extracted based on a variety of aperture 
sizes on the purpose of aperture correction. The final magnitudes 
with aperture correction are calibrated with photometric standard 
stars, which are HD 110621 for /-band (/ = 8.91), HD 188112 
for //A"-band (H = 10.78, K = 10.89), and HD 161743 for L- 
band (L = 7.61). The astrometry is adjusted by referencing the 
NACO point sources with the detections of the SINFONI data 
whose astrometry was calibrated using the 2MASS catalog. With 
this procedure, the relative astrometric difference between the 
NACO data and SINFONI data is better than 071. 

The TIMMI2 (Thermal Infrared Multimode Instrument; 
iReimann et al.l iT998l f mid-IR imaging was carried out in 2003 
July at the ESO 3.6 m telescope at La Silla, Chile. The ob¬ 
servations covered the VI, A'10.4, and Q 1 bands with ,Ln of 
8.7 pm, 10.38 /mi, and 17.72/mi, respectively, and all had FOV 
of 55" x 38" with pixel scale of 072. All data are limited by 
diffraction with a FWHM (Full Width Half Maximum) of 077. 

The M 17 UC1 -IRS5 region was imaged with VISIR (VLT 
Imager and Spectrometer for mid-InfRared'. ^Lasase et akll2004l) 
in 2006 May through the Sic filter (4 cl i = 11.85 pm). The obser¬ 
vation procedures are described in N+07. The image is of good 


quality (FWHM ~ 0732) limited by diffraction, with pixel scale 
of 07127. The astrometry of the TIMMI2 and VISIR imaging 
data were calibrated on the basis of the NACO data. The astro- 
metric accuracy is better than 071 throughout the FOV centered 
on the M 17 UC1 -IRS5 region. 

The spectroscopy of the /V-band silicate absorption feature 
was performed with TIMMI2 at the ESO 3.6 m telescope at La 
Silla, Chile, within the same observation run of UC1 (N+07). 
The seeing was 077; the slit width was 172. 

4. Results 

4. 1. High angular Resolution Near-IR imaging 

The M 17 UC1-IRS5 region outlines a typical interface be¬ 
tween an H ii region and a PDR, as characterized by a prominent 
IF seen in our previous , lower resolution near-IR images (e.g. 
iHoffmeister et al.l 120 08. CZ+12). Our high-resolution images 
unveil many fine structures toward this region which were not 
revealed before (see Fig. |T}. Part of the results had been pub¬ 
lished for UC1 (N+07). Interestingly, the two luminous infrared 
objects (M 17 IRS5, M 17 B273) are resolved to have more than 
one component in the high-resolution near-IR images. If they ac¬ 
tually are multiple systems, one needs to review these objects be¬ 
cause they were treated as single high-mass YSO candidates, and 
thus the primary’s brightness would have been overestimated. 

The JHKL magnitudes of the objects inside the white box 
are listed in Table [2] The brightest component of the IRS5 sys¬ 
tem is treated as the primary star, namely IRS5A, and followed 
by four other fainter companions IRS5B, IRS5C, IRS5D, and 
IRS5E. Similarly, B273A is chosen as the primary star of the 
binary B273. 

4.1.1. Point sources with IR excess 

Rg. m shows the HKL color-color diagram of the objects from 
Table[2] Objects located to the right of the reddening vectors are 
suggested to show infrared excess which traces the circumstel- 
lar material around YSOs. The level of infrared excesssomehow 
relates to the amount of circumstellar material dLada & Adams! 
Il992h . which in turn reflects the evolutionary stage of YSO. All 
the objects with accurate L-band magnitudes show infrared ex¬ 
cess, indicating that they are YSOs. IRS5A, B273A and B273B 
all show infrared excess resembling to classical T Tauri stars 
(CTTS). IRS5B seems to have a larger infrared excess, imply¬ 
ing a younger evolutionary stage. IRS5C is the reddest source 
apart from UC1, although its infrared excess is comparable to 
that of IRS5B. Without proper spectral type classification, we 
can only establish a crude mass sequence for these point sources 
according to their dereddened magnitudes. However, it is impos¬ 
sible to extract exact extinction values for these sources solely 
based on their locations in Fig. [2] The close positions in the 
HKL color-color diagram suggest that all sources except IRS5C 
have similar reddening. The extinction of IRS5C, however, is 
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Fig. 1 . Three-color image of the M 17 SW PDR taken by NACO at three near-IR broadband filters (blue: H ; green: K; red: L), with the area of the 
SINFONI integral field spectroscopy outlined in the white box. Individual point sources discussed in thi s paper are labeled. Three types of maser 
are marked according to th eir coordinates: 22-GHz wa ter masers (filled triangles, [Caswe ill 9981: 1 Johnson et alll998h : Class II methanol masers at 
6.66 GHz (filled diamonds. ICaswell et al.|[l995l . l2000h : OH masers at 1.67 GHz (thick plus. ICasweilTT998fk 

Table 2. NACO JHKL photometry of the point sources in M 17 UC1-IRS5 region 



Source 

ID 

R.A. 

(J2000) 

DEC 

(J2000) 

J 

(mag) 

H 

(mag) 

K 

(mag) 

L 

(mag) 

M 17 B273 

A 

18 20 25.08 

-16 11 34.0 

14.03 + 0.01 

12.19 + 0.01 

10.71+0.01 

8.89 + 0.01 

B 

18 20 25.14 

-16 11 34.6 

20.11+0.33 

17.18 + 0.05 

15.26 ± 0.04 

13.25 + 0.09 

M 17 UC1 

18 20 24.84 

-16 11 35.1 

- 

- 

13.61 + 0.02 

6.20 ± 0.03 

M 17 IRS 5 

A 

18 20 24.60 

-16 11 39.4 

13.34 + 0.01 

11.38 + 0.02 

9.81+0.01 

8.16 + 0.03 

B 

18 20 24.54 

-16 11 39.0 

19.07 + 0.22 

15.09 + 0.02 

13.81+0.02 

11.56 + 0.03 

C 

18 20 24.54 

-16 1141.2 

- 

19.13 + 0.13 

14.65 ± 0.03 

9.69 + 0.01 

D 

18 20 24.63 

-16 11 40.7 

- 

17.97 + 0.12 

16.32 + 0.07 

15.27 ± 0.40 

E 

18 20 24.69 

-16 11 40.5 

- 

19.69 + 0.32 

17.75 + 0.18 

- 


Notes. Coordinates of point sources adopted here are a compromise between the astrometry of NACO and SINFONI data, whose differential has 
been checked less than 071 throughout the FOV. 


about three times higher. From the //-band magnitudes, we sug¬ 
gest that IRS5A is most massive source, followed by B273A and 
IRS5C, IRS5E has the lowest mass while the other sources have 
masses in between. 

The area enclosed by the white outline harbors several types 
of maser, such as water (H 2 O) maser at 22 GHz, hydroxyl (OH) 


maser at 1.665 GHz and Class II methanol maser at 6.66 GHz. 
Such a variety of maser is c ommonly found toward regions 
of ma ssive star formation (e.g. iFontani et al.ll2010L iBreen et aD 
120101) . Particularly, the Class II methanol maser is believed to 
be radiatively pumped, and typically coincides in position with 
hot molecular cores, UCH 11 regions, OH masers and near-IR 
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Fig. 2. NACO K - L versus H - K color-color diagram. The blue curve 
denotes the intrinsic colors of main-sequence stars dDucati e t al.|[200 1[): 
the cyan line corresponds to the locus of T Tauri stars lMeveretal.1 
119971) . Re ddening vectors (dashed l ines) with slope of E{H - K)/E(K - 
L) = 1.2 (Hoffmeister et al. 2008) are drawn. Their lengths equal to 
visual extinction of 90 mag. 

sources JFontani et al.ll2010l) . Note that FFO, OH, and 6.66-GHz 
methanol masers are all detected in the close vicinity of UC1, 
which coincides with the expectation of a massive forming star. 
Moreover, a 6.66-GHz methanol maser is detected in close prox¬ 
imity to IRS5A. The H 2 O maser at 22 GHz is closer to IRS5C 
than to IRS5A. If the 22-GHz HiO maser is associated with 
IRS5C, it might be a high-mass YSO too. The other possibility 
is that the 22-GHz H 2 O maser is associated with IRS5A, coex¬ 
isting with the 6.66-GHz methanol maser. 

4.1.2. Diffuse emission 

A bright rim crossing between the B273 binary is clearly seen 
in Fig. Q] This rim is a portion of the IF in M 17 seen in larger 
view (e.g. CZ+12) that represents the boundary between the H 11 
region (to the northeast) and M 17 SW (to the southwest). We 
note an elongated, bar-like emission feature 2" west to the IF. 
This emission bar is nearly parallel to the IF, and shows similar 
brightness and color with the IF. However, the emission bar has a 
smooth boundary while the IF is very sharp. Hence, the emission 
bar might not have the same physical origin as the IF. 

The near-IR polarization studies toward M 17 (CZ+12) re¬ 
vealed a bar-like feature enhanced in AT-band polarized light ly¬ 
ing west to the IF. Fig.0shows the comparison between AT-band 
polarization and NACO /K image. The polarization pattern cen¬ 
tered on IRS5A had been discussed in CZ+12. The high polar¬ 
ization degree of UC1 is consistent with its two reflection lobes 
(N+07). At the position of the emission bar seen in NACO /K im¬ 
age, a concentration of polarization vectors ~ 10% - 15% marks 
the location of an IRN. From the pattern of these polarization 
vectors, UC1 is the best candidate for the illuminating source of 
the IRN. In fact, the coincidence the two bar-like features seen 
in different manners indicates a mass concentration of gas and 
dust associated with UC1 in a bar-like shape. 



Fig. 3. IRSF/A - polarization vectors (red lines; CZ+12) overlaid on 
NACO/A" image (grey scale). 


4.2. Mid-IR imaging 

Previous mid-IR imaging showed round morphologies both for 
UC1 and IRS5 with angular resolution 1 - 2". Interestingly, the 
new VISIR image at 11.85pm with angular resolution ~ 0'.'3 
reveals substructures seen in Fig. 0 for both objects. The IF is 
still visible at this wavelength, but much fainter than in the K- 
band. 

UC1 shows an elongated feature toward northwest, which 
agrees well with its circumstellar disk orientated at P.A. (from 
north to east) « 146° (shown in yellow dashes in Fig. 0 . More¬ 
over, the two parts separated by the circumstellar disk are asym¬ 
metric, with the northeastern part more extended. This asym¬ 
metric structure can be explained by an inclined circumstellar 
disk with angle » 30° with respect to the line of sight (here¬ 
after LOS), which is proposed by N+07 to reproduce the K -band 
asymmetric structure which possesses a much brighter south¬ 
western lobe. Besides the mid-IR emission in the close vicinity 
of UC1, a structure extending to north is more worthy to note. 
Interestingly, it is spatially coinciding with the IRN discussed in 
Sect. 14.1.21 and fully covers the area of substantial K -band po¬ 
larization. Conversely, the bar-like feature seen in AT-band covers 
only a portion of the IRN. 

The high-sensitivity VISIR image shows two lobes separated 
by a dark lane orientated at P.A. ~ 30° for IRS5A. Unlike UC1, 
the Af-band polarization level is very small in the bulk area of 
IRS5A. And the A"-band polarization found in east and south of 
IRS5A has no mid-IR counterparts; only the mid-IR feature in 
northeast of IRS5A coincides with A"-band polarization. 

Besides the two bright objects mentioned above, B273A and 
IRS5C are also visible in the VISIR 11.85 yum image; IRS5C is 
even brighter than B273A at this wavelength. 

In addition, the mid-IR images at other wavelengths taken 
by TIMMI2 show only the two bright mid-IR objects, UC1 and 
IRS5A, because of the lower angular resolution and sensitivity. 
Thus, these mid-IR images taken by TIMMI2 are not shown 
here, but only the flux densities of UC1 and IRS5A are listed 
in Table 0 
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Table 3. Mid-IR flux densities of the YSOs in the M 17 UC1-IRS5 region 


Object 

8.7 pm a 
(mJy) 

10.38 pm a 
(mJy) 

11.85 pm 
(mJy) 

17.72 pm" 
(mJy) 

20.6 pm c 
(mJy) 

M17 B273A 

- 

- 

46 ± 22* 

- 

- 

M17UC1 

18700 ±1300 

7300 ±1000 

31300+ 1100" 

146700 ± 29700 

128500 ± 6800 

M17IRS5A 

3200 + 1600 

6800 ±1500 

9700 ±1100" 

130000 + 31000 

103900 ± 5500 

M17IRS5C 

- 

- 

256 + 25* 

- 

- 


Notes. (n) 270 aperture. (i> 074 aperture. (c) 372 aperture from K+02. 
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Fig. 4. VISIR 11.85 pm image of the M 17 UCHRS5 region. IRSF/K 
polarization vectors (blue lines; CZ+12) are overplotted as well. The 
disk position angle is denoted by the yellow dashes with RA of 156° 
(N+07). 
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Fig. 5. Three-color image created from the emission (line+continuum) 
of three near-IR lines (red: Bry 2.167 pm, green: H 2 1 -0 S(l) 2.122 pm, 
blue: He 1 2.059 pm). The three squares (I, II, and III) mark selected 
positions representing Hu region, ionization front, and PDR. A strip 
line with tickmarks is also drawn perpendicular to the IF of M 17 SW. 
Offset at each tickmark is -1, 0, +2, +4, +6, +8, +10, +12 in arcsec 
from northeast to southwest, respectively. 



4.3. SINFONI near-IR spectroscopy of diffuse content 

The SINFONI observations are centered on IRS5, and are ex¬ 
tended northeast to cover UC1 and B273. Fig.[5]shows the three- 
color composite of three lines. This figure reveals both the stellar 
content and diffuse emission in that area. For each object visible 
in Fig. [5J a SINFONI H + K spectrum is available. However 
a spectral classification could only be obtained for IRS5A and 
B273A (Sect. ELl, while the S/N of the other spectra was not 
sufficient. Bry and 2.059 pm Hei 2 1 P-2 1 S line emission are de¬ 
tected throughout the FOV. The location of the H 11 region and 
the IF is traced by the very strong Bry emission at the top-left 
corner. On the other hand, the PDR is characterized by 2.122 pm 
H 2 1-0S(1) emission. We note Ho emission surrounding IRS5A 
and Ho emission coinciding with the bar-like IRN illuminated by 
UC1. 

4.3.1. Nebular emission lines 

The spectra in the three selected regions are shown in Fig. [6j and 
the observed lines are listed in Table [4] The hydrogen recombi¬ 
nation lines such as Brackett series and Pa a are visible in all re¬ 
gions. Besides the prominent atomic hydrogen lines, five atomic 
helium lines are also observed; among them 2.059 pm Hei is 


the strongest one. Molecular hydrogen emission lines longward 
2 pm are also observed in region II and III. An emission fea¬ 
ture at 2.287 pm is generally classified a_s_Hi v = 3 - 2S(2) 
line or unidentified (UID) line (e.g. lLumsden et aI1l200ll) . How¬ 
ever, the Hi/UID line might have been mismatched for the H 2 
v = 3 - 2 S(2) line, because H 2 emission is not spatially co¬ 
incident with this H 2 /UID line in planetary nebulae (e.g. NGC 
7027, lOkumura et al. 120041) . The most potential carrier of this 
UID line is [Se iv] at 2.287 pm, which is typically detected in 
highly excited ISM suc h as planetary nebulae and UCH 11 regions 
with hot O-type stars dBlum & McGregon l2008, and references 
therein). The ionizing sources (04 binary) of M 17 H 11 region are 
very likely to produce the 2.287 pm Se iv line in M 17 SW. The 
2.287 pm UID line keeps roughly the same strength throughout 
the FOV of Fig. [3 indicating that this line traces the ionized gas 
region. Therefore, we attribute the 2.287 pm line to [Se iv]. 

4.3.2. H i/H 2 transition zone of M 17 SW 

The coexistence of Bry and H 2 1 — 0S(1) in region II and III 
implies that region II and III maybe part of the H 1 /H 2 transition 
zone of M 17 SW. In order to investigate the scale size of this 
H 1 /H 2 transition zone, Fig. [3 shows line strength variations of 
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A (//. m) 

Fig. 6. Ambient spectra of the three selected positions marked by squares in Fig. [5] The telluric feature in the range 1.8 - 2.1 pm is not corrected, 
because its strength is in proportion to the nebular continuum, which is much weaker than the emission lines. Thus the near-IR emission lines are 
little affected by the telluric feature. 


Bry, 2.059gin He i and H 2 1—0 S( 1) along the strip line denoted 
in Fig. [5] Peaks of Bry and 2.059/rm Hei both occur around 
IF. Specifically, the 2.059/rm Hei peak occurs ~ O'.'l closer to 
the H ii region than the Bry peak, which defines the IF. Molec¬ 
ular hydrogen starts to appear at offset ~ +0"4, closer to the 
PDR. Along the strip line, the line strengths of Bry and 2.059 /rm 
He i show almost identical variations, which both experience fast 
growth until the IF, rapid drop (offset < +0'/5), shallow decline 
(offset between +075 - +2"), and roughly constant baseline (off¬ 
set longwards +2"). The H 2 1 — 0S(l) line strength has another 
type of variation, which shows four peaks along the strip line. 
The coexistence of atomic gas (H i. He i) and molecular gas (H 2 ) 
in the range +0.4"-+11.7" and the trend further into the PDR in¬ 
dicates that the H i/H 2 transition zone of M 17 SW almost starts 
from the IF and extends toward the cloud core with a projected 
scale size more than 12" (0.12 pc). However, one has to be cau¬ 
tious because the emission from the HII region and the PDR may 
partly overlap along the LOS if the configuration is not strictly 
seen edge-on; the H n region emission lines originate from the 
surface of the PDR, while PDR emission lines form inside the 
PDR. In such a case, the co-existence of emission lines like 
in Fig. |7| will not tell us the true story about Hi/H 2 transition 
zone. Comparison with a strictly edge-on PDR will help to clar¬ 
ify whether this scale size of H i/H 2 transition zone in M 17 SW 
is reliable or not. 


The Orion Bar is a dense PDR like M 17 SW, but nearly 
edge-on (lAllers et alJ l2005[ and references therein). We mea¬ 
sured the scale size of Hi/H 2 transition zone in the Orion Bar 
according to Fig. 3 of lHavashi et alJ (119851) to be around 20", 
which is consistent with the distributions of H 2 rotational lines 
dAllers et ali 12005 ). Due to the edge-on configuration of the 
Orion Bar, the scale size measured above is identical to the real 
scale zie of the Hi/H 2 transition zone, i.e. 0.04 4 pc at the dis¬ 
tance of Orion (450 pc: iHoogerwerf et~akl l 200Q ;). However, the 
Hi/H 2 transition zone in M 17 SW is at least three times Orion 
Bar in size, which can not be simply explained by the differ¬ 
ential properties between the two PDRs. We speculate that the 
presence of the three emission lines in Fig.|7]is simply a result of 
geometric projection. Unlike the Orion Bar, M 17 SW is indeed 
inclined from the LOS with a substantial angle. This result con¬ 
firms the conclusion of a recent study about the LOS structure 
of M 17 SW based on Hi rotational emission (ISheffer & Wolfir e 
120131 ) . 

4.3.3. H 2 excitation 

As seen in Table [3 six H 2 emission lines are detected. The two 
most prominent H 2 lines are H 2 1 — 0Q(1) at 2.407 pm and 
H 2 1 - 0 S( 1) at 2.122 pm. H 2 emission lines in PDR are gen¬ 
erally thought to have two physical mechanisms - fluorescent 
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Table 4. Observed near-IR lines toward the M 17 UC1-IRS5 region 


A (pm) 

IDs 

Transition 

Eupper(K) fl 

F±6F(10 

-15 -2 -1 

erg cm s 

1 arcsec z ) 

Obs. 




region I 

region II 

region III 

1.5349 

HI 

18-4 


4.3 + 2.0 

4.9 ±2.1 

1.1 ±0.4 

1.5447 

HI 

17-4 


4.2 ± 1.5 

4.8 ± 1.6 

1.7 ±0.3 

1.5566 

HI 

16-4 


3.7 ± 1.2 

2.6 ± 1.0 

1.8 ± 0.6 

1.5710 

HI 

15-4 


3.3 ± 1.1 

4.9 ± 2.2 

2.1 ± 0.2 

1.5891 

HI 

14-4 


5.4 ± 0.6 

5.7 ±0.7 

2.1 ± 0.6 

1.6120 

HI 

13-4 


8.4+ 1.9 

7.6 ±2.6 

3.9 ± 1.2 

1.6417 

HI 

12-4 


6.8 ± 2.1 

5.6 ±2.6 

4.5 ± 0.4 

1.6817 

HI 

11-4 


11.6 + 0.8 

11.4 + 0.5 

7.0 ± 0.5 

1.7013 

He I 

4 3 D-3 3 P 


7.3 ± 1.6 

7.4+ 1.6 

3.3 ±0.3 

1.7372 

HI 

10-4 


17.6+ 1.1 

16.8 ± 1.0 

9.9 ± 0.6 

1.8185 

HI 

9-4 


23.3 ± 1.3 

21.8 ± 1.8 

13.0+1.3 

1.8696 

He I 

4 3 F-3 3 D 


7.9 ±0.9 

7.8 ±0.8 

5.7 ± 0.7 

1.8762 

HI 

4-3 


183.9 ± 1.0 

147.8 ± 1.1 

74.8 + 0.7 

1.9457 

HI 

8-4 


31.7+ 1.0 

32.1 ± 1.5 

19.5 + 0.8 

2.0344 

h 2 

v =1-0 S(2) 

7584 

- 

0.6 ± 0.4 

1.8 ±0.4 

2.0592 

He I 

2 1 P-2 I S 


61.2 + 0.6 

60.8 + 0.6 

28.4 + 0.3 

2.1132 

He I 

4 i. 3 _ 3 i. 3 p 


4.2 ± 0.8 

4.0 ± 0.3 

2.4 ± 0.4 

2.1224 

h 2 

v =1-0 S(l) 

6956 

- 

3.3 ±0.4 

8.5 ± 0.2 

2.1620 

He I 

7 3 F-4 3 D 


3.1 ±0.4 

3.0 ±0.5 

1.5 ±0.4 

2.1667 

HI 

7-4 


99.4 ± 0.7 

92.8 + 0.4 

57.3 + 0.3 

2.2239 

h 2 

V =1-0 S(0) 

6471 

0.0 ± 0.7 

1.5 ±0.5 

3.2 ± 0.6 

2.2873 

[Seiv] 

2 P3/2— 2 Pl /2 


1.2 ±0.4 

1.4 ±0.4 

1.5 ±0.5 

2.4073 

h 2 

v =1-0 Q( 1) 

6149 

- 

7.2 ± 1.0 

15.7 + 0.5 

2.4141 

h 2 

v =1-0 Q(2) 

6471 

- 

- 

4.8 ± 0.4 

2.4244 

h 2 

v =1-0 Q(3) 

6956 

- 

4.1 ± 1.2 

6.8 ± 0.7 


Notes. Lines detected in the avera ged spectra of three selected regions. The flux of each line is averaged over each region. 
Excitation energy adopted from lDabrowskii dl984i) . 



Offset (arcsec) 

Fig. 7. Line intensity variations of three representative lines - Bry, 
2.059 pm He i, and H 2 1 - 0 S( 1) - along the strip line which is per¬ 
pendicular to the IF of M 17 SW (outlined in Fig- [5j- The zero-point 
is defined at 18: 20: 25.136, -16: 11: 34.21 (J2000), at the maximum 
strength of the IF; negative offsets point toward the Hu region, while 
positive offsets point toward the PDR. 


excitation by FUV photons and thermal excitation in shock 
fronts. Practically, near-IR H 2 fluorescent spectra have been ob¬ 
served for a variety of classical PDRs such as associated with 
NGC 2023, Orion Bar, and the northern bar in M 17 ( iGatlev etal . 


l!987l:lHavashi et al.l 1985L [Tanaka et al.ll989l) . Meanwhile, a va¬ 
riety of H 2 spectra have been found to be produced in the shock 
fronts associated with the jets/outflows found in PDRs (e.g. 

Niirnberger et al.lf2007t iMartm-Hernandez et al.li2008t: iBik et al.l 

201 Oh . Although these two mechanisms both act occasionally for 
H 2 emission found in PDRs, they can be distinguished since the 
corresponding H 2 line ratios are different. 

A common way of characterizing the H 2 emission is to eval¬ 
uate the gas temperature in the framework of a ro-vibrational di¬ 
agram, which is a plot of the observed column density against the 
energy of the upper level. The column density, Nj, of the upper 
level of a given transition can be calculated from the measured 
line intensity, /, of the corresponding H 2 line via the following 
formula: 


Nj = 


47r V 

Ajhc ' 


where Aj (the rest wavelength) and Aj (the Einstein A- 
coefficient) are taken from iTurner et al.l d 1977b . If collisional de¬ 
excitation plays dominant role, the H 2 molecule will be in LTE 
and the energy level obeys the Boltzmann distribution. In such 
distribution, the relative column densities of any two excitation 
levels can be expressed in terms of excitation temperature T ex : 


Ni gi -( E t - Ej ) 

— = — exp - 

Nj gj kT ex 


where g, is the degeneracy, Ej is the excitation energy taken 
from iDabrowskil d 1984b . and k is the Boltzmann constant. If the 
gas is thermalized at a single temperature, the plot of the log¬ 
arithm of Nj to gj ratio against the level energy will reveal a 
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Fig. 8. Column densities of Hi 1-0 S(l), 1-0 Q(l) and 1-0 Q(2) (quadrants I to III) as well as the Hi excitation temperature map (quadrant IV). 
In quadrants I to III, IRS5A.C and UCl’s positions are marked by asterisks. In quadrant IV, the ro-vibrational diagram based on the H 2 lines 
(Hi 1-0S(1), 1-0S(0), 1—0Q(l), 1-0Q(2). and 1 —0Q(3)) averaged over all collected pixels is plotted, as well as the histogram density of Hi 
excitation temperature. 


straight line, whose slope provides the reciprocal of the excita¬ 
tion temperature. 

We implied this method to evaluate the excitation tempera¬ 
ture for the detected Hi emission lines, assuming the same dust 
attenuation for all lines. Fig.[8]shows the extinction-uncorrected 
column densities of Hi 1-0S(1), Hi 1-0Q(1), and H 2 1-0 Q(2). 
We excluded the Hi 1 —0 S(2) line in plotting the ro-vibrational 


diagram, because this line is too weak to be significantly detected 
in most FOV. Except H 2 1-0S(1) and 1-0Q(l), the remaining 
three Hi lines are weak, although stronger than 1-0 S(2). To in¬ 
crease the S /N ratios for these weak Hi lines, we degraded the 
angular resolution of the SINFONI data by rebinning the data 
with 3x3 array. The resulting S /N ratios are increased by a factor 
of - 3. Note that rebinning is applied merely when the Hi emis- 
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sion lines are analyzed. Thus, only Fig. [ 8 ] and Fig. [TO] are plot¬ 
ted in the rebinned pixel scale. For H 2 1-0 S(l) and 1-0Q(l), a 
threshold of S/TV > 3 is used; for the remaining three H 2 lines, 
S/TV > 2 is applied. The lower-right panel in Fig. [ 8 ] presents 
the map of T ex for H 2 emission, and the histogram of excita¬ 
tion temperatures, as well as the ro-vibirational diagram based 
on the column densities averaged throughout the collected pix¬ 
els. The ro-vibrational diagram shows that all five FF lines lie 
along a straight line whose slope corresponds to an excitation 
temperature ~ 600 K. We note that the determined 7 K for each 
rebinned pixel might have large errors because the adopted FF 
lines just cover a very narrow range of excitation temperature 
(6100 - 7000 K). Nevertheless, the statistical T ex of the FF gas 
is more meaningful than a specific value for each rebinned px- 
iel. With the histogram of T ex for all rebinned pixels, we found 
a peak around 575 K. A single excitation temperature for the FF 
emission indicates that the energy levels are thermally excited. 
Other discriminators for the two excitation mechanisms of FF 
emission such as FF 1 —0S( 1 )/2— 1 S(l) ratio and FF ortho-to- 
para ratio both indicate thermal excitation. 

For fluorescence emission, each vibrational level has Njigj 
lying along a separate “branch”, and the rotational population 
within each level can be approximated by a thermal distribution. 
If there are transitions from several levels, a curved line in the ro- 
vibrational diagram would therefore provide evidence for non- 
LTE process dMartm-Hernandez et ali [2008!) . The FF lines at 
higher levels, e.g. v = 3-2, can be useful to distinguish between 
the two mechanisms. In the case of the M 17 UC1 -IRS5 region, 
the merely available v = 1-0 transitions cannot be used to 
discriminate between shock and fluorescent excitation, because 
they will lie along a straight line in the ro-vibrational diagram 
for both mechanisms if uncertainties are considered. Moreover, 
in dense PDRs the lower FF levels will be thermalized as like in 
shock fronts. This degeneracy between dense PDRs and shock 
fronts prevents us to exactly characterize the FF emission from 
the ro-vibrational diagram. This approach is more meaningful 
for the regions with moderate/strong FF emission at all transi¬ 
tions. Indeed, it is curious to know whether weak FF emission 
area possessing the same excitation mechanism with the stronger 
FF emission area. 

In an H ii region, another attempt to qualitatively character¬ 
ize the FF emission comes from the comparison of molecular 
and atomic hydrogen emission. The molecular-to-atomic line ra¬ 
tio, particularly the ratio FF 1-0 S(l)/Bry, is < 1 in active galax¬ 
ies (starburst, Seyfert, and ultra-luminous infrared galaxies) and 
star-forming regions (lHatch et al.ll2005l and refere nce therein). 
Conv e rsely, this ratio is > 1 in ou tflow regions (lHatch et alJ 
12005b iMartm-Hernandez et alJ 1 20081) . All regions showing FF 
1-0 S(l) emission yield FF 1-0S(1)/Bry « 1. The extremely 
strong hydrogen atomic emission observed here indicates very 
intense incident stellar radiation field, which makes the FUV 
pumping as the most likely excitation mechanism for the molec¬ 
ular hydrogen emission observed in our study. Further indirect 
evidence also points out that FUV fluorescence is the most plau¬ 
sible excitation mechanism. For instance, FF 1—0S(l) and FF 
1-0Q(1) emission seen in the top two panels in Fig. [ 8 ] shows 
some filament structures parallel to the IF. This FF emission can 
not be characterized by the aforementioned approaches, since 
just emission of 1-0 S(l) and 1-0 Q(l) are available. One may 
argue, however, that the FF emission lying parallel to IF can be 
produced in shock waves driven by the IF. An identical case is 
the Orion Bar PDR which shows FF emission behind itself with 
a ratio of H? 1 —0 S( 1 )/2-1 S(l) ~ 3 r esembling sho c k exci tation 
(lHavashi et alJ 1 19851) . Nevertheless, IBurton et alJ (|T990) pro- 



18 1, 20 m 25.2 s 25.0 s 24.8 s 24.6 s 24.4 s 

Right Ascension (J2000) 


Fig. 9. FF 1—0 S(l) line flux (red contours) overlaid on the IRSF/Tf po¬ 
larization degree image (grey scale; CZ+12). FF 1—0S(l) contour lev¬ 
els start from 2 x 10 -15 erg cnT 2 s _1 arcsecT 2 with four intervals, each of 
2 x 10~ 15 erg cnr 2 s” 1 arcsec~ 2 . Green contours correspond to the VISIR 
11.85 pm emission. 


posed a high-density PDR model to explain the observed line ra¬ 
tio because the typical shock speed, < 3 km s , of an expanding 
H ii region driving into molecular gas is too low to significantly 
excite low-v FF transitions. Considering the common properties 
between M 17 SW and Orion Bar, we tend to neglect shock exci¬ 
tation for the IF-parallel FF emission. Complementally, in Fig. [9] 
the IRN illuminated by IRS5A and by UC1 (the one discussed 
in Sect. 14.1.2b spatially coincides with FF emission. FF fluores¬ 
cent emission has been found for some large-scale IRN associ- 
ated with star-form i ng regions such as N GC 2023 and Orion Bar 
(iGatlev et al .11 19871 IBurton et alJll99C)h . Although the IRN here 
are associated with luminous IR stars, we speculate that this co¬ 
incidence is an indirect evidence in favor of FUV fluorescence 
excitation. 

In Fig. [TO] we show the Fla 1—0 Q(l)/1—0S(1) ratio map, 
which covers most of the Ha emitting regions. The Ha 
1 —0 Q( 1 )/l —0 S( 1) ratio intrinsically varies with the physi¬ 
cal conditions of Ha gas. iLuhman et al.l (i!998l) modeled Ha 
line ratios detected in Orion Bar and Orion S including Ha 
1 —0 Q( 1)/1 —0 S( 1) ratio for several set of PDR models and ther¬ 
mal excitation by shock front as well. From their modeling, we 
find distinct Ha 1 —0 Q( 1 )/l —0 S(l) ratios between PDR models 
and shocked thermal excitation (see Table 1 in iLuhman et al.l 
H998h . The H 2 1-0Q(1)/1-0S(1) ratios of the H 2 emission in 
M 17 SW are mostly in the range 1.0 - 3.0. The bulk of the H 2 
emission regions have a H 2 1 —0 Q( 1 )/l —0 S( 1) ratio ~ 1.8, which 
is between the ratios of two PDR models of mh = 10 6 cm~ 3 , Go = 
10 4 , T 0 = 500 K, and n H = 10 6 cm~ 3 , G 0 = 10 5 , T 0 = 1000 K, re¬ 
spectively. In contrast, the H 2 1 —0 Q( 1 )/l —0 S( 1) ratio predicted 
by the shocked thermal model is 0.7, much lower than the ratios 
observed here. The extinction-corrected H 2 1 —0 Q( 1)/1 —0 S( 1) 
ratio will be smaller than the observed one, however, the de¬ 
crease caused by dereddenning is not big enough to lower this 
ratio to a thermal value. 

We note quite uniform distributions of line ratios such as H 2 
1—0 S(l)/2—1 S(l), H 2 1-0S(1)/Bry, and H 2 1-0 Q(1)/1-0S(1) 
across all H 2 emitting areas in the M 17 UC1-IRS5 region, 
which indicates a uniform mechanism for the entire H 2 emis- 
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Right Ascension (J2000) 

Fig. 10. Map of H 2 1—0 Q(l)/1—0 S(l) ratio, and corresponding his¬ 
togram density at lower left corner. The three asterisks are as same as in 

Fig. 0 


sion. All the above properties of H 2 emission can be well ex¬ 
plained in the scenario that the H 2 molecules inside M 17 SW 
are initially pumped to vibrational states by FUV fluorescence, 
and then are deexcitated from high-v levels to low-v levels due 
to the high-frequency collision inside the high-density PDR. Due 
to this process, H 2 emission is mostly in the form of low-v vibra¬ 
tional states and pure rotatio nal states. Besides the v = 1 - 0 H 2 
emission lines reported here. ISheffer & Wolfire! (1201 3l) reported 
mid-IR pure-rotation H 2 emission at v = 0-0 states toward M 17 
SW. 


4.4. SINFONI near-IR spectroscopy of point sources 

The SINFONI H+K spectra of the YSOs in M 17 UC1 -IRS5 re¬ 
gion, are crucial to understand their fundamental properties such 
as effective temperature and luminosity class via analyzing the 
characteristic spectral features (see Fig. ITil. Objects located in¬ 
side an HII region or a PDR require special care due to the con¬ 
tamination by nebular emission. The nebular contamination is 
evaluated in an annulus closely encircling the stellar contours. 
The spectra of IRS 5A and B273A have sufficient S /N to main¬ 
tain a reliable classification, while the other sources are too faint. 

The most prominent lines in Fig. [IT] are the hydrogen Brack¬ 
ett absorption lines seen in IRS5A, the CO 2-0 bandhead ab¬ 
sorption for B273A, and the Hi Paa and Bry emission for UC1. 
The Hi emission lines of UC1 are reminiscent of the common 
emission features of YSOs, which are due to the accretion flows 
falling onto the YSOs. However, we note from Fig. 0 that the 
emission lines mostly form in the two reflection lobes of UC1, 
and very little part directly comes from UC1 which is mostly 
obscured by its edge-on disk. The Hi emission lines seen in 
UCl’s spectrum arise most likely in the associated HCHII re¬ 
gion. The shape of H 1 emission region oriented perpendicularly 
to the edge-on circumstellar disk of UC1, is suggested to trace 
the relevant HCHII region expanding preferentially along the po¬ 
lar direction because of the lower-density gas in polar regions. 


The normalized spectra of IRS 5A and 273A are shown 
in Fig. [12] and [13] respectively. IRS5A’s spectrum displays 
strong hydrogen Brackett absorption lines. The sole charac¬ 
teristic line in the K -band, Bry, indicates a temperature class 
later than B3, since its effective temperature is not high enough 
for the growth of 2.113 pm Hei line which appears for early- 
B stars ( Hanson et al.lll996 T The Bry equivalent width (EW) is 
quantified to disentangle the ambiguity between mid-B (kB4- 
B7) and late-B/early-A (kB8-A3); the latter shows Bry EW 


greater than 8 A while the former shows Bry EW between 4 A 
and 8 A (Hanson et al. 1 9% ). The Bry EW of IRS5A is mea¬ 
sured to be 4.8 A, constraining its temperature class to be kB4- 
B7. How does the K -band temperature class_link to the optical 
spectral class? In the study by lHanson et al.l (119961) . six refer¬ 
ence stars with known optical spectral type are assigned to a K- 
band spectral type of kB4-B7. Four out of the six stars have op¬ 
tical spectral classes of B3-B6, with three dwarfs and one giant. 
The remaining two are optically classified late-B/early-A super¬ 
giants. jHanson|etjiU (1996) stressed that the spectral types solely 
based on K -band spectrum are not sensitive to surface gravity 
and show ambiguity between mid-B stars and late-B/early-A su¬ 
pergiants. 

The //-band spectrum provides more spectral type indica¬ 
tors. The presence of strong Bril line and weak 1.701 /imHei 
line (EW < 0.2 A) both point out a spectral class later than early- 
B (lHanson et al.l 19981) . Moreover, the Brl 1 EW of IRS5A is 
determined as 4.1 A, suggesting mid-B spectral class. With the 
known spectral type range of IRS5A, it is possible to distin¬ 
guish dwarfs and s upergi ants based on the Brl 1/1.701 pm He 1 
ratio. lHanson et al.l (1 19981) found for early/mid-B stars that this 
ratio is consistently larger for dwarfs than for supergiants of the 
same spectral class; e.g. this ratio is about 3.0 in average for B4- 
B7 supergiants. In the case of IRS5A, the large ratio of Brl 1 to 
1.701 /imHe 1 (> 20) indicates a luminosity class close to dwarf. 


Some Hei lines (e.g. 1.701yum, 2.113/im) are crucial diag- 
nostic lines to separate early-B sta rs from later type stars (e.g. 
lHanson et al.l l2005: Bik et al.|l2005l) . Moreover, the difficulty in 
distinguishing dwarfs and giants solely based on near-IR spectral 
features also prevent us to better constrain the luminosity class 
of IRS5A. According to the H + K spectral characteristics, we 
suggest a spectral type of B3-B7 V/III for IRS5A. 

B273A’s H -band spectrum does not show any significant 
characteristic line. Its //-band spectrum only shows 2.3 pm 
CO 2-0 bandhead absorption. The CO bandhead absorp- 
tion longwards of 2.29 pm is typical of late-type stars (e.g. 
[Wallace & Hinklelll997l) . whose outer atmospheric layers have 
the proper temperatures (~ 1000 - 3000 K) to produce such 
features. Moreover, low-mass Class II/III YSOs (temperatures 
identical to those of late-type stars) might also show CO band- 
head in absorption despite of the circumstellar material be¬ 
cause the circumstellar veiling is less than that of Class I 
YSOs and consequently overwhelmed by the photospheric fea¬ 
ture ( iCasali & Eiroafl996l) . 


The CO bandhead absorption is found to be tightly related 
with the effective temperature and luminosity class for giants and 
supergiants. The EW of CO bandhead increases linearly with 
temperature declining for giants and supergiants, respectively, 
and the later has higher EW than the former of the same effec¬ 
tive temperature (e.g. lGonzalez-Fernandez & Negueruelal2012i) . 
The CO bandhead of B273A coincides in depth and width with 
that of the reference star, a G4 giant (see Fig. o. Indeed, a su¬ 
pergiant needs to have earlier temperature class to match the 
EW of a giant’s CO bandhead. We measured the EW of CO 
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Fig. 11. Flux-calibrated SINFONI H + K spectra of the point sources in M 17UC1-IRS5 region. All natural spectral features are marked, e.g. Hi 
emission/absorption lines and CO 2-0 bandheads in absorption. 
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Fig. 12. Normalized SINFONI H + K spectrum of IRS5A (black curve). For comparison, t he spectra of several reference stars are displayed; 
they are adopted from the literature as follows: B2 V (HD 36166), AOV (HD 146624) from lHanson et all ll2005t) : B2 V/IV (HD 19374), B3III 
(HD 209008), B8 V (HD 169990) from lHanson et all (t 1996T) . All reference spectra are converted to a resolution identical to SINFONI's resolving 
power. 


bandhead for B273A to be « 3.5 A, which suggesting a spec¬ 
tral type of G4/G5 for a giant and G2/G3 for a su pergiant when 
following Fig. 2 in iGonzalez-Fernandez & Negueruelai d2012l) . 
Alternatively we checked the EW of CO bandhead for the 
main sequence stars catalogued bv l Wallace & Hinkle (19971) and 
found that B273A’s CO bandhead strength is between that of a 
G8 V (HR 4496) and a K2 V (HR 1084) star. We could not further 
constrain the spectral type of B273A without other characteristic 
lines at shorter wavelengths. 

A fit example is CEN 34, a mid-G supergiant located along 
the LOS toward M 17, is classified to be background with re¬ 
spect to M 17 with age between 50-100 Myr ( IChen et aT1l2() 13l) . 
Upon the consideration of B273A’s apparent brightness at short 
wavelengths, the object should not be in the foreground of M 17. 
Thus we can rule out the possibility of an early-K dwarf which 
is only visible at shorter distance. On the other hand, the solu¬ 
tion of early-G supergiant puts B273A in the backside of M 17. 
However, the large amount of molecular gas along the LOS to¬ 
ward the M 17 UC1 -IRS5 region can almost obscure any back¬ 
ground object. Therefore, we keep the solution of mid-G giant 
as the spectral type of B273A. 

4.5. Extinction and Luminosity 

In order to place IRS5A and B273A in the Hertzsprung-Russel 
diagram (HRD), one has to determine their luminosity. In prin¬ 
ciple, the already constrained spectral type of the source would 
give some hints to derive the extinction toward the source, be¬ 
cause the temperature-dependent intrinsic color is reddened by a 
certain amount of extinction to the observed color. For instance, 
a general manner of determining the extinction for B stars is ap¬ 
plying the short wavelength colors (e.g. V-R vs I-R) rather than 
the near-IR colors (i.e., / - H vs H - K ), on consideration that 
the former ones are more sensitive to the temperatures of early- 
type stars. In the cases of YSOs, the short wavelength colors are 
the only reliable way to derive the extinction value because the 
near-IR colors are usually affected by the strong near-IR excess 
emission which is emitted by the hot circumstellar material. 

The absence of hydrogen emission lines in the spectra of 
IRS5A and B273A suggests that there is no substantial accretion 
activity. We speculate that the J - H color is not severely af- 



A (/im) 


Fig. 13. Normalized SINFONI spectrum of B273A (black) in the wave¬ 
length range 2.27 - 2.42 pm, covering CO 2-0 handhead absorption. 
For comparison, the reference spectrum of mid-G giant HIP 53316 (SpT 
G4 Ilh lWallace & Hinkle! 19971) is shown in red. The reference spectrum 
(R ~ 3000) has been smoothed to the spectral resolution of SINFONI 
spectrum. 


fected by dust thermal emission due to the absence of hot accre¬ 
tion tunnel flows. Following the extinction law at near-IR from 
ICardelli et ak (1 19 891) and the s pecific value of Ry = 3.9 toward 
M17 from iHoffmeister et akl (120081) . the Ay values of IRS5A 
and B273A are estimated from their apparent / - H colors on 
the basis of the intrinsic J - H colors at certain temperatures 
dPecaut & Mamaiek 1201 3l) . The reliability of the obtained ex¬ 
tinction is tested by reddening a blackbody of the proper temper¬ 
ature with the same Ay value. Note that the lCardelli et al.l f 1989 ) 
near-IR extinction law is invalid in the 5 - 20 pm range, hence 
the extinction law used by iRobitaille et al.l (2 0071 ) that works in 
the 1 - 20pm range is employed when reddening the SEDs. In 
the 1 — 3 pm range, the ex tinction law used bv IRobitaille et all 
d2007l) is equivalent to the ICardelh et al.l (1 19891) extinction law 
when Ry = 3.1. The reddened blackbody is shown in blue in 
Fig. [14] The blue blackbody curve matches the bluer part of 
the stellar SED for both objects. On the other hand, IRS5A dis- 
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plays a deep 9.1 pm silicate absorption feature. Using the em¬ 
pirical r elation for silicate optical depth Ay = (18 + 1 ) 79.7 
(Whittetj[2003l) . IRS5A’s A v is estimated as ~ 41 mag from 
Tgj ~ 2.3. Note that Ay of IRS5A is even higher than that 
(~ 40 mag) of UC1 (N+07), while the latter is redder than the 
former. Such a contradiction even could not be explained by the 
difference in spectral ty pe between them. Alter natively, if we 
adopt Ay/rg. 7 = 9-19 dRoche & Aitkenlll985l) . Ay of IRS5A 
is hence between 21-44 mag. The Ay based on Tgj has a large 
error which originates from the various Ay /Tgj ratios. The Ay 
value of IRS5A calculated from its J - H color agrees with the 
upper limit from optical data. IRS5A and B273A’s Ay values de¬ 
rived by their J-H colors are adopted in this paper (see Tabled. 

The visual extinction Ay evaluated above is then converted 
to Ay at /-band following the relation Ay ~ 0.27 Ay when 
Ry - 3.9. Further, the /-band absolute magnitude ( Mj ) can 
be estimated for both IRS5A and B273A through distance 
m odulus. For IRS5A, the va lue of Mj is -3 mag. Referring 
to IPecaut & Mamaiek] (12013l) . the temperature of IRS5A (~ 
15000 K) corresponds to a /-band bolometric correction factor 
(BCy) of -1.6 mag. Hence, the bolometric magnitude, Mboi, of 
IRS5A is estimated to be -4.6 mag, corresponding to a bolo¬ 
metric luminosity ~ 5500 L Q . Here the solar M^i is adopted as 
+4.75 mag. Similarly, we have Mj » -l.Omag, and BCy » 1.25 
for B273A (L e ff ~ 5200 K). Mbd ~ 0.25 of B273A is equal to a 
bolometric luminosity of 63 L (:) . 

The flux measurements from near- to mid-IR wavelengths 
provide an observational approach to evaluate the luminosity. 
With the new flux measurements in the 1-20 pm range pre¬ 
sented in this paper, we are able to update the IR luminosity of 
IRS5A, and derive the IR luminosity for IRS5C and B273A. As 
shown in Fig. [14] we employ the simple blackbody curve to fit 
the infrared flux of the three objects. The significant absorption 
dip around 10 pm in IRS5A’s SED corresponds to the 9.7 pm 
silicate absorption which is confirmed by the TIMMI2 spectrum. 
Hence the flux measurements at wavelengths around 9.7 pm are 
not used in fitting the SED of IRS5A. However, for the SEDs of 
IRS5C and B273A, because the 11.85 pm flux is the only avail¬ 
able measurement in mid-IR, we still use this flux measurement 
in fitting. The 9.7 pm silicate absorption should be taken into ac¬ 
count when dereddening the SED. Be aware of that the 11.85 pm 
flux is affected by the 9.7 pm silicate absorption feature, thus 
the blackbody fitting will overestimate the color temperatures of 
IRS5C and B273A. 

The SED of IRS5A can be well fitted by two blackbody 
components with color temperatures of 1204 K and 264 K. The 
cooler component indicates a large amount of dust grains along 
the LOS of IRS5A. The color temperatures of IRS5C and B273A 
are 496 K and 1218 K, respectively. The observed IR luminosity 
in the 1 - 20pm range, Lir, of these objects are summarized in 
the third column of Table 0 

IRS5A’s Ljr is very close to the values in N+01 and K+02, 
which reported values of 3000 L e and 2600 L (: , for IRS5, respec¬ 
tively. In their work, the aperture size is larger than the one 
used in this paper, and IRS5C is not resolved in their larger 
beam infrared images. A smaller luminosity derived for IRS5A 
in this paper is reasonable. On the other hand, the observed Lir 
of IRS5A and B273A are much lower than their bolometric lu¬ 
minosity Lboi derived above. If the observed flux measurements 
are extrapolated to longer and shorter wavelengths, and are cor¬ 
rected for the extinction, the L ho | of IRS5A and B273A can be 
constrained observationally. 

The extinction-corrected Lir is obtained for IRS5A and 
B273A by using the proper extinction law mentioned above. Fur- 




Fig. 14. SEDs constructed based on the flux measurements summarized 
in Tableland Table [3] (filled circles), and the M-band point by ISAAC 
(filled square; ITT()ffnieisterll20()8l) . For IRS5A, the TIMMI2 spectrum at 
around 9.7 pm is denoted by the solid curve in the top panel. The tem¬ 
perature of the blackbody (in blue) corresponds to the proper spectral 
type according to lCoxI 12000 ). 


thermore, the Lboi of IRS5A and B273A are derived by extrap¬ 
olating the best-fitting blackbody to longer and shorter wave¬ 
lengths, e.g. in the range 0.02 - 200 pm. The obtained Lboi is 
listed in the fourth column of Table 0 

The Lboi of IRS5A and B273A derived from SED fitting 
are compared with the values obtained from the temperature- 
dependent bolometric correction values. For IRS5A, the former 
value is significantly larger than the latter one. This difference 
is partly caused by the extinction used in dereddening the SED 
of IRS5A. IRS5A shows evidence of circumstellar material, the 
extinction caused by the circumstellar material is included when 
dereddening the SEDs. However, the circumstellar extinction 
does not attenuate the Lboi, but just redistributes the SED be¬ 
cause the circumstellar material absorbs the blue photons and re¬ 
emits at longer wavelengths if assuming an ideal blackbody. In 
the case of IRS5A, the circumstellar extinction might be a small 
portion of the total extinction, upon the consideration that the 
foreground extinction toward the southwestern bar is su ggested 
to be on the order of 10 mag or even higher (e.g. 'Pellegrini et al.l 
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Table 5. Infrared luminosity of the YSOs in the M 17 UC1-IRS5 re¬ 
gion 


Object 

Ay 

(mag) 

Lir* 

(Lq) 

^bol 

(L Q ) 

M17B273A 

13 

10 (44) 

71 

M 17 UCl fa 

40 

2.5 x 10 3 (1.1 x 10 4 ) 

3.0 x 10 4 

M 17 IRS5A 

18 

2.0 x 10 3 (3.0 x 10 3 ) 

8.6 x 10 3 

M 17 IRS5C 

- 

13 

- 


Notes. 

Values in brackets are extinction-corrected Lir. 
(6) Values adopted from N+07. 


120071: ISheffer & Wolfirdl2013h . We recalculated the L ho | with a 
smaller Ay value by 1 mag difference, and got smaller L h() i of 
5.8 x 10 3 L 0 for IRS5A, which agrees much better with the value 
of 5.5 x 10 3 L e . Therefore, we take the Lboi values in Table[5]as 
the upper limits, and the extinction-corrected Lir as the lower 
limits in the following discussions. 

5. Discussion 

5.1. Multiplicity 

With the high angular resolution infrared images toward M 17 
UC1 - IRS5 region, the companions of the luminous infrared ob¬ 
jects (IRS5, B273) are resolved for the first time. According to 
the HKL color-color diagram, the two components of the B273 
system both show L-band excess emission. Thus, we suggest that 
B273A and B273B form a young binary system with a projected 
separation of about 2200 AU. Three of the five components of 
IRS5 show L-band excess emission. In light of IRS5D’s position 
in Fig. [2] we also regard it as a member of IRS5 system taking 
into account the photometric errors, and suggest IRS5 might be 
a young quadruple system. 

5.2. Structures of the circumstellar material 

The large mid-IR excess emission of IRS5A and UC1 previ¬ 
ously reported by N+01 and K+02 indicate circumstellar ma¬ 
terial for both objects. In addition, from a near-IR polarization 
view (CZ+12), the high-level polarization individually associ¬ 
ated with IRS5A and UC1 also traces dust grains physically 
linked with the two objects. Specifically, the almost edge-on cir¬ 
cumstellar disk of UC1 found by N+07 produces the high-level 
/L-band polarization vectors parallel to the disk. The polariza¬ 
tion pattern of the IRN illuminated by IRS5A, however, implies 
a quite different orientation than UC1. E.g. the inclination angle 
of IRS5A’s IRN estimated in CZ+12 is ~ 50° from the plane of 
the sky. 

Fig. [9] presents perfect geometric match between the Hi 
emission and polarized light, which explains the physical car¬ 
riers of the polarized light. The FF emission residing just be¬ 
side UC1 is exactly aligned with the elongated feature extending 
to northwest, which is suggested to coincide with the dust ther¬ 
mal emission originated in the circumstellar disk of UC1 (see 
Sect. 14. 2I >. We speculate that such FF emission is related with the 
circumstellar disk of UC1, and actually traces the warm FF gas 
in UCl’s disk. The dust grains in UCl’s disk scatter the incident 
starlight, and produce the polarized light seen at near-IR. The 
FF emission filament related to UCl’s disk is one of the three FF 
emission filaments that are parallel to each other, and parallel to 
the IF. The FF filament at the midway between the IF and UC1 


is spatially coinciding with the bar-like emission feature seen in 
the NACO/Zr image (see Sect. 14. 1 .21 and Fig. [3}. Moreover, the 
area covered by this bar-like emission in turn is also occupied by 
the south-north mid-IR extended emission originated from UC1. 
This south-north extended emission seen in mid-IR and the co¬ 
ordinated IRN seen in near-IR polarized light together suggest a 
south-north outflow from UC1. The bar-like emission seems to 
arise at the northeastern side of this outflow because this side is 
facing to the IF. The IF northeast of the outflow suppresses the 
outflowing material on the northeastern side, and consequently 
forms a sheet of gas with higher temperature than the rest of the 
outflow. The temperature difference between the sheet at out¬ 
flow’s northeastern side and the rest of the outflow explains well 
the different morphology of the outflow seen between near-IR 
and mid-IR wavelengths. Moreover, the coordinated emitting FU 
gas at the northeastern side of the outflow is suggested to be 
excited at a temperature of ~ 600 K (see Sect. 14.3.31) . while 
the rest of the outflow is suggested to have thermal tempera¬ 
ture around 100 K heated by an early-B outflow driving object 
dLongmore et al.ll201 It) . The north-south outflow of UC1 is most 
likely the blue-shifted lobe, with the red-shifted lobe invisible 
because of the severe attenuation. 

The coordinated emitting FF gas surrounding IRS5A and the 
IRN illuminated by the same object indicates a distribution of 
gas and dust embracing IRS5A in the north with an arc-like 
shape. The circular shape of the IRN illuminated by IRS5A re¬ 
minds us that the southern half would be associated with FF gas 
either, although the southern half is not fully covered by the SIN- 
FONI data in its small FOV. The temperature of this material dis¬ 
tribution is ~ 600 K, similar to the associated emitting FF gas. 
With the high-resolution VISIR 11.85 pm image, we expect a 
tight relation between the diffuse mid-IR emission and the IRN, 
just like what has been observed for UC1. However, the diffuse 
mid-IR emission is only observed at the midway between IRS5 A 
and UC1, where merely a portion of the IRN is spatially associ¬ 
ated (see Sect. 143 . Indeed, this mismatch can be the result of 
two separate material distributions that are just along the same 
LOS; the dust temperature responsible for the mid-IR emission 
of IRS5 A is modeled to be around 264 K, while the dust temper¬ 
ature linked with the IRN and FF emission is fitted to be around 
575 K. The two lobes in the core region of IRS5A stack with 
the inner circle of the arc-like material along the LOS. The dif¬ 
ferent morphology of IRS5A seen at near-IR and mid-IR rise 
the question if the mid-IR emission comes from the area close 
to IRS5A, or from its envelope or outflow. The dark lane sep¬ 
arating the two mid-IR lobes can be somehow explained by a 
nearly edge-on cold disk that absorbs the mid-IR emission from 
backside. The proposed cold disk should have impacts at near-IR 
bands, such as a silhouette dividing the object into two lobes as 
well, very similar to UC1. Nevertheless, IRS5A’s near-IR im¬ 
ages do not show any sign of silhouette disk. The extinction 
of IRS5A is another strong issue against the assumption of a 
nearly edge-on cold disk, because the circumstellar extinction is 
only a small portion of the total extinction of IRS5A. We pro¬ 
pose that the mid-IR emission actually arises from the outflow 
whose axis is almost along the LOS of IRS5A. The two lobes 
seen in mid-IR correspond to the outflow walls at two sides. In 
view of the severe attenuation at near-IR bands and at around 
10/rm ( Chapman et al.l l2009i> . the Hj emission filament is the 
blue-shifted lobe. Meanwhile the outflowing material traced by 
the mid-IR emission might be blue-shifted, because of the atten¬ 
uation issue too. 

We note from Fig. |8]that IRS5C is closely associated with H 2 
1-0S(1) emission. Indeed, IRS5C is embraced by Hi 1 —0S( 1) 
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Fig. 15. Positions of the YSOs in the HRD with well constrained temperatures and luminosities. Red rectangles outline the parameter space of 
IRS5A (the upper one) and B273A (the lower one), and the red asterisk at upper left marks UC1. For IRS5A and B273A, the conversion from 
spectral types to effective temperatures and bolometric luminosities is referred to Sect. 1431 while UC1 is assumed to be a ZAMS B0.5 star in light 
of the HCHII region (N+07 and references therein). Left panel: positions of the YSOs compared with the model of massive star formation via disk 
accretion (H+10). The symbols along the evolutionary loci at various accretion rates and the ZAMS locus denote the mass at each position. Right 
panel: positions of the YSOs placed against the PMS evolutionary tracks of different masses (CA+12). Isochrones of five ages are drawn as the 
dashed lines in different colors. 


emission which shows an elliptical shape orientated at P.A. 
~ 65°. The bulk of H 2 emission is shifted northeast to IRS5C. 
As described in Sect. 14.1.11 IRS5C is coinciding with 22-GHz 
HjO maser which is located 075 northwest to IRS5C. The 
H 2 1-0S(1) emission associated with IRS5C resembles that 
found for some Herbig Ae/Be stars in their circumstellar disks, 
where molecular H 2 gas is heated by the Ae/Be stars them¬ 
selves dCarmona et al.ll201lh . The strong excess emission at L- 
band (this work) and mid-IR indicates the presence of circum¬ 
stellar material surrounding IRS5C. Hence it is plausible that 
H 2 gas among the warm circumstellar material thermally emits 
H 2 1-0 S(l) emission detected here. However, unlike the Hi 
l-OS(l) emission detected in the circumstellar disks of HD 
97048 and HD 100546 which originates at radii on the order of 
~ 10 1 - 10 2 AU, the size of the H 2 1—0 S(l) emission associated 
with IRS5C is about 1600 AU in diameter. H 2 emission occur¬ 
ring at large distances requests massive stars other than Herbig 
Ae/Be stars to maintain sufficient radiation field. Nevertheless, 
the other possibility is that the H 2 emission is externally excited 
by the FUV photons from the ionizing OB stars in M 17. In light 
of the excitation mechanism of the H 2 emission associated with 
other objects, the latter explanation is more plausible and reli¬ 
able. To explain the location of the H 2 emission, we propose a 
nearly edge-on flared disk with radius of 1600AU to approxi¬ 
mately match the shape of H 2 emission. In addition, the 22-GHz 
H 2 0 maser resides 075 away along the axis of the circumstellar 
disk of IRS5C. In some cases, the H 2 0 masers can be excited 
in the shock fronts where the outflowing gas from high-mass 
YSO strongly interacts with the surrounding molecular gas (e.g. 
iTorrelles et all2014l) . Considering the position of the H 2 0 maser 
with respect to the proposed disk, and the disk’s size, we inter¬ 
pret that IRS5C is a deeply embedded high-mass protostar with 
a nearly edge-on flared disk. IRS5C does not show radio free- 


free emission at 1.3 cm (iJohnson et alJl99~8i) . which implies that 
IRS5C is still at an early stage. 

5.3. Evolutionary Stages of YSOs 

With the known temperature and luminosity for IRS5A and 
B273A, the locations of these two objects in the HRD will re¬ 
veal their evolutionary stages when compared with numerical 
models for star formation. The model of H+10 specifically con¬ 
siders the evolution of a massive protostar via disk accretion, 
while that of Iciaretl d2012h models the PMS evolution purely 
tuned by the gravitational contraction without mass feeding, i.e. 
non-accretion model. 

Both models indicate that IRS5A is still experiencing PMS 
stage. Using the accretion model H+10 to estimate IRS5A’s 
mass and age is problematic, due to the constant accretion rate 
and fixed initial condition used in H+10. If referring to the non¬ 
accretion PMS model, IRS5 A’s mass is currently about 9 M e and 
will not change much when IRS5A is evolving toward ZAMS. 
And a non-accretion 0.1 Myr isochrone passes by IRS5A’s posi¬ 
tion in the HRD. Note that non-accretion model simplifies the 
evolution of protostars, whose accretion histories, initial radii 
and thermal efficiencies can strongly alter their positions in the 
HRD. Nevertheless, very young age estimates (< 1 Myr) for stars 
hotter than 3500 K based on non-accretion model are reliable 
dHosokawa et al.ll201 ll) . Therefore, we suggest that the age esti¬ 
mate for IRS5A is reliable. 

However, the circumstellar material encircling IRS5A that 
is found to be the bipolar molecular outflow driven by itself, 
and the lack of accretion tracers both doubt that IRS5A would 
continue accretion to gain mass. The peculiarity of this object 
is that it currently does not show any tracer of accretion funnel 
flow, which imply either periodic accretion or feedback (radia¬ 
tion pressure, outward mass flow) halted accretion. Both mecha- 
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nisms are thought to have impact on the final mass of high-mass 
protostar, hence are crucial to understand the formation of high- 
mass protostar. Otherwise, the detection limit of the SINFONI 
data could also prevent resolving weak hydrogen emission lines 
(see the discussion in Sect. l5.4b . Follow-up investigations for the 
gas kinematic in the close vicinity of IRS5A will contribute to 
determine the gas motion and consequently clarify the evolu¬ 
tion trend via constraining the accretion status. Such follow-up 
studies need high spatial resolution and high spectral resolution 
observations, which can be implemented by ALMA (Atacama 
Large Millimeter/submillimeter Array). 

B273A is located in the right-lower area of HRD. Model 
H+10 interprets mass gaining for B273A, however, it seems to 
lack hydrogen recombination lines. The excess emission at L- 
band and mid-IR both indicate the presence of circumstellar ma¬ 
terial, e.g. a circumstellar disk. Hence, we propose phases of pe¬ 
riodic accretion or halted accretion for B273A, while the emis¬ 
sion line strength below the detection limit could be plausible 
too. In the other case, the location of B273A in the right panel 
of Fig. [15] suggests a mass ~ 4 M 0 and an age ~4x 10 5 yrs. As 
a similar case to IRS5A, we believe that B273A’s age estimate 
based on the non-accretion PMS model is reliable, and B273A 
is still at the first-half stage of its entire evolution life toward 
ZAMS. 


5.4. Feedback on accretion 


One may argue that the spectra of IRS5A and B273A have a 
low resolution (R ~ 1500) and not high enough S/N, which 
leads to the ’fake’ non-detection in case of weak Bry emission. 
Based on the spectral resolution (FWHM ~ 3 pixels) and con¬ 
tinuum’s noise level, we estimated that the upper limit of Bry 
emission flux is ~ 2.2 x 10" 15 ergs cm" 2 s" 1 for IRS5A, and 
7.9 x 10 16 ergs cm" 2 s 1 for B273A; with the extinction values 
of this two objects shown in Table [5] the extinction-corrected 
Bry line flux, / Bry , is ~ 1.3 x 10" 14 ergscm~ 2 s"' for IRS5A, 
and ~ 2.8 x 10 15 ergs cm -2 s 1 for B273A. Further the lumi¬ 
nosity of Bry emission line was computed as L Bly = 4 n d 1 / Bry , 
where d is the distance of M 17. The Bry line luminosity is found 
to be tightly correlated with the accretion luminosity, L acc , for 
a stat i stical sample of low - mass YSOs (e.g. Mendieutfa et alJ 


1201 it iRigliaco et alJ I20T2L lAlcala et alJ 2014 ). For high-mass 
YSOs, the correlation between L Bly and L acc is not well estab¬ 
lished due to the lack of a statistics meaningful sample. As an 
approximation, we used the log L Bly - log L acc relation reported 
in I Alcala et ali (120141) to estimate the L acc , which is ~ 2.2 L (: . for 
IRS5A, and ~ 0.38 L 0 for B273A. From the positions of this two 
objects in the HRD, we estimated a radius of ~ 10 R 0 for both 
of them. Together with the masses of the two objects mentioned 
above, the formula 


Hacc 


1.25 


L a cc R* 
G M„ 


gives the upper limit of M acc , which is ~ 1.1x10 7 M 0 yr 1 for 
IRS5A, and ~ 3.8 x 10" 8 M 0 yr" 1 for B273A. 

High-mass stars can arrive at ZAMS stage still during their 
accretion phase. However, the enormous radiation pressure and 
ionizing photons will act on the accretion flow and dissipate 
the circumstellar material. Thus, the timescale of the accretion 
phase for a high-mass protostar will strongly depend on the stel¬ 
lar mass, which in turn is constrained by accretion rate. Such 
timescale is far to be known for high-mass protostars. UC1 ’s disk 
suggests that this B0.5 ZAMS star might be still in its accretion 


phase. A plausible evidence supporting accretion comes from the 
Bry emission at the position of UCl’s southwest lobe (Fig. [5j. 
Nevertheless, Bry emission can be produced in the HCHII re¬ 
gion associated with UC1. For an HCHII region, the dynamics 
and morphology of the ionized gas are determined by the ratio of 
the ionization radius, and the gravitational bound radius, R/,. 
defined by G M*/ 2 c 2 , where M, is the stellar mass and c s is the 
sound speed dKetoll2007l) . With M„ = 15 M e and c s = lOkms" 1 
(at the temperature of 10 4 K for ionized gas), R/, is 66 AU. Since 
the Bry emission must come from the ionized gas of the HCHII 
region, the Bry line map in Fig. [5] shows the morphology of the 
HCHII region around UC1. The distance from the obscured cen¬ 
troid of UC1 to the southwestern Bry emission knot is taken as 
the value of R h which is ~ 0725, or ~ 500 AU. Thus, for UC1, 
R, » R/,. corresponding to a stage of HCHII region that the 
ionized gas moves outward to form an outflow ( Kclo 2007 ). In 
this picture, the accretion is confined to a narrow range of angles 
close to the mid-plane of the disk and is close to termination. 
Without the observational evidence of the gas dynamics in the 
inner most region of UC1, the accretion status of UC1 is contro¬ 
versial. 

Due to the large inclination angle of IRS5A from the plane 
of the sky, near-IR emission lines produced in the accretion flow 
(if exists) can be seen without the severe attenuation by the disk. 
In fact, IRS5A does not show any emission line indicating accre¬ 
tion. We have proposed two reasons for the absence of accretion 
indicators in Sect. 15.31 If in the case of non-accretion phase, the 
evolutionary stage of IRS5A would imply that accretion could 
even be halted before high-mass protostar arriving at ZAMS. In 
view of the lower temperature at the early stage of high-mass 
protostar, radiation feedback is not as strong as for a ZAMS star, 
thus outflow activities are additionally required to be responsible 
for the reversed accretion. In a large picture, the expanding H n 
region of M 17 might also contribute to the dissipation of accre¬ 
tion tunnel flows. The deduced timescale of accretion phase is 
thus ~ 1 x 10 5 yrs, on the same order of high-mass protostar’s 
age. On the other hand, the probable periodic accretion process 
of IRS5A would suggest a timescale of accretion phase much 
longer. The determination of the gas dynamics in the inner most 
area of IRS5A will disentangle this confusing puzzle. 

Interestingly, the absence of hydrogen recombination lines 
in B273A’s spectrum implies a terminated-accretion phase. An 
intermediate-mass YSO with an age ~ 4 x 10 s yrs like B273A 
is most likely to show accretion activity, because YSOs at later 
evolutionary stages such as intermediate-mass T Tauri stars 
(1.5 - 4M 0 , 1 - lOMyr) were observed to show Bry e missi on 
even in low-resolution (R ~ 800) near-IR spectra (ICalvet et al.1 
|2004| ). This discrepancy between the observation and the predi¬ 
cation for B273A indicates additional feedback to dissipate the 
accretion funnel flows. The exact spatial coincident between 
B273A and the IF implies that the expanding Hu region driven 
by the ionizing photons is the plausible mechanism of external 
feedback. 

From the current masses and ages of these two ob¬ 
jects, IRS5A and B273A require mass accretion rate of 9 x 
10" 5 M 0 yr" 1 and 1 x 10" 5 M 0 yr" 1 , respectively. Even in the pos¬ 
sibility of weak Bry, we still note a dramatically drop of M for 
both objects. This could be also the result of feedback from stel¬ 
lar activities and expanding H n region. For the two low-mass 
YSOs, B273B and IRS5B, we could similarly estimate that the 
upper limit of their mass accretion rates is definitely lower than 
thatofB273A,i.e. ~ 3.8 xlO _8 M 0 yr _1 . The typical mass accre¬ 
tion rate of low-mass YSOs is on the ord er of 1 x 10 8 M 0 yr" 1 
(e.g. iFang et al.ll2(H3L lAlcala et al.ll2014l) . which is comparable 
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to the detection limit of B273B and IRS5B’s spectra. The lack of 
Bry emission line in the spectra of low-mass YSOs (e.g., B273B, 
IRS5B) is reminiscent of non-accretion phase, however, we can¬ 
not rule out the possibility that their accretion indicators could be 
revealed by higher resolution spectroscopic observation. In the 
case of non-accretion phase for these low-mass YSOs, it seems 
to remind us that such external radiation feedback might be im¬ 
portant in the disk clearing process for the YSOs in massive star 
forming regions. 

5.5. Properties of M 17 Sl/V 

Previous studies of M 17 SW speculated that a dense PDR 
with high-density clumps and low-den sity gas to match 
the observed propert ies of M 17 SW (iMeixner e t al.l ll 992t 
IPerez-Beau nuits et akll2010l 120121) . However, these studies are 
based on far-IR and radio line observations, and thus cannot re¬ 
solve structures smaller than several arcseconds. The near-IR 
H 2 emission detected in the M 17 UC1 -IRS5 region with ex¬ 
tremely high spatial resolution (100 mas) set a entirely new pic¬ 
ture for the warmest part of M 17 SW. The column densities 
of several Hi v = 1-0 ro-vibrational lines show that H 2 gas 
mostly concentrates in filaments in the range 0.02 - 0.05 pc. 
The H 2 line ratios, as been discussed in Sect. 14.3.31 refer to a 
gas density of_> 10 5 crrT 3 if they are compared to PDR mod¬ 
els dBurton et al.l 11990 ). The high-density nature of M 17 SW 
is also evidenced by the molecular gas persisting just west of 
the IF (iPerez-Beaunuits et al.ll2012l) . This can be interpreted as 
the result of self-shielding of molecules in high-density PDR 
which can move the transition regions for Hi/H 2 close to the 
surface of the cloud dBurton et al.il 1 99(1) . Complementally, the 
H 2 gas temperature characterized by the ro-vibrational diagram 
is in the range 500 - 750K. In contrast. Hi observations to¬ 
ward M 17 specifically showed that Hi number density toward 
the M 17 UC1-IRS5 region is in the range 3.7 x 10 3 - 1.5 x 
10 4 cm 3 , which corresponds to Hi spin temperature between 
50 K and 200 K in many directions toward M 17 dBrogan et alJ 
119991 ). Therefore, the H 2 filaments presented here confirm for 
the first time that dense clumps with number density two or¬ 
ders of magnitude higher and temperatures higher than ambient 
atomic gas can exist down to 1 0~ 2 pc scale inside a dense PDR 
like M 17 SW, which was previously suggested being clumpy 
dMeixner et al.lfl992h . 

6. Summary and conclusions 

In this paper we presented diffraction-limited near- to mid-IR 
images and SINFONI integral field spectroscopy at H + K to¬ 
ward the M 17 UC1 -IRS5 region, which shows emission lines 
of Hei. Hi, and H 2 . Our diffraction-limited data reveal new fine 
structures of this region. 

1. This work complements the work by CZ+12 from a view 
of much higher angular resolution. The IRN identified in 
CZ+12 illuminated by IRS5A and by UC1 are confirmed 
to trace the molecular outflow driven by the two objects, 
respectively. Combining the SINFONI H 2 line emission map 
and the mid-IR VISIR image, a blue-shifted outflow lobe is 
proposed for IRS5A, with two types of tracers in the form of 
H 2 emission filament and mid-IR emission. The molecular 
outflow UC1 is running south-north with its blue-shifted 
lobe merely visible in forms of near-IR polarized light and 
mid-IR emission. The northeastern side of UCl’s outflow 
is suppressed by the shock fronts driven by the expanding 


Hn region. Thus a sheet of warm, dense gas forms at the 
northeastern edge of UCl’s molecular outflow, and appears 
in forms of the bar-like emission seen in K -band and the H 2 
emission filament which are both parallel to the ionization 
front. 

2. The uniform line ratios of all H 2 emitting areas indicate the 
same excitation mechanism; the H 2 molecules are initially 
pumped by FUV photons, and are repopulated by the colli- 
sional deexcitation in a dense PDR. The H 2 gas excitation 
temperature is estimated to be around 575 K based on the 
ro-vibrational diagram of the detected H 2 v = 1 lines. The 
Hi line ratios are used to probe the properties of the PDR. 
Comparisons to PDR models suggest clumpy PDR structure 
with at least two components: warm, dense molecular 
clumps with «h > 10 5 cm 3 and T ~ 575 K and atomic 
gas with ~ 3.7x 10 3 - 1.5x 10 4 cm' 3 and T ~ 50-200K. 

3. IRS5 might be a young quadruple system containing IRS5A, 

IRS5B, IRS5C, and IRS5D. Its primary star IRS5A is clas¬ 
sified as a B3-B7 V/1II star with a bolometric luminosity of 
3.0 - 8.6 x 10 3 confirmed to be a high-mass protostar with 

mass ~ 9 M q and age ~ 1 x 10 5 yrs, while the three lower 
mass companions are much less constrained. Particularly, 
IRS5C might be a deeply embedded high-mass protostellar 
object with a dusty disk. The spectral type of B273A is 
assigned to G4/5 III, suggesting an intermediate-mass YSO 
of ~ 4 M q and an age ~4x 10 5 yrs when compared with the 
non-accretion PMS model. B273A may have a lower mass 
companion B273B of possibly similar age. 

4. The Bry emission of UC1 is ambiguous for tracing the accre¬ 
tion because the emission can arise from the HCHII region as 
well; UC1 might have terminated accretion from a concern 
on the ratio of R , » R/,. The absence of Bry emission in the 
SINFONI spectra of the other YSOs of various masses im¬ 
plies terminated accretion for them. Even in the case of weak 
Bry emission that cannot be resolved by the SINFONI data, 
the steep drop in the accretion rates of IRS5A and B273A 
along the protostellar evolution suggests processes of dissi¬ 
pating accretion funnel flows. The object IRS5A is unique 
because it is during a stage when its radiation feedback is 
not high enough to severely destroy accretion funnel flows 
or even to terminate the accretion; therefore, the expand¬ 
ing H ii region and its bipolar outflow both are suggested to 
contribute simultaneously. For the other YSOs, the expand¬ 
ing Hu region could be the major mechanism on dissipat¬ 
ing/terminating accretion. 
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